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THEORY OF RARE B DECAYS 



A. Ali 

Deutsches Elektronen- Synchrotron DESY, Notkestrasse 85, D-22603 Hamburg, ERG 

We discuss some selected topics in rare B decays in the context of the standard 
model and compare theoretical estimates with available data. Salient features of 
the perturbative-QCD and power corrections in the decay rate for B — > Xs + 7 
are reviewed and this framework is used to determine the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix element \Vts\, yielding \Vts\ = 0.033 ± 0.007 from the 
present measurements of the electromagnetic penguins. We review estimates of 
the ratio Rk' = r{B — > K* + 'y)/r{B — > Xs +7) in a number of theoretical 
models, which give a consistent account of this quantity. Issues bearing on the 
photon energy spectrum in _B — > Xs + 7 are also discussed. The CKM-suppressed 
decays B — > + 7, B^ — > + 7, and — > (p^jOi) + 7 are reviewed with 
particular emphasis on the long-distance contributions in the exclusive decays. 
The impending interest in these decays in determining the parameters of the CKM 
matrix is emphasized. Finally, the semileptonic decays B Xs£'^i~ are also 
discussed in the context of the SM. 



1 B{B Xs+j) and B ^ K* +j in the Standard Model and exper- 
iment and determination of \Vts\ 

1.1 Experimental status 

Electromagnetic penguins were first sighted on the B territorji in 1993 by the 
CLEG collaboration through the exclusive decay B K* + 7EI. This feat was 
followed by the daunting measurement of the inclusive decay B Xg + 7 in 
1994 by the same collaborationla. The present CLEG measurements can be 
summarized asB: 

B{B^Xs+j) = (2.32 ± 0.57 ± 0.35) X 10"'*, 
B{B~^K*+j) = (4.2±0.8±0.6) X 10~^ (1) 

which yield an exclusive-to-inclusive ratio: 

i....^||^f±^Hl8.1±6.8)%. (2) 

Very recently, the inclusive radiative decay has also beeiL reported by the 
ALEPH collaboration with a (preliminary) branching ratio □: 

B{Hb X,+j)^ (3.29 ± 0.71 ± 0.68) x IQ--*. (3) 

Since the ALEPH measurement is done at the peak in the process — > 
bb ^ Hi) + X ^ (Xg -f 7) -f A, the branching ratio in (0) involves a different 



1 



weighted average of the various _B-mesons and A;, baryons produced in 
decays (hence the symbol Ji\, ) than the corresponding one given in (|^), which 
has been measured in the decay T(4S') ,B^B^ . Theoreticahy, the 

inclusive radiative decay widths for the various beauty hadrons are expected to 
be nearly equal. Despite this, their branching ratios are not all equal reflecting 
the differences in the respective total decay rates (equivalently lifetimes). 

In the context of SM, the principal interest in the decay rates in eqs. 
and (^) lies in that they determine the ratio of the CKM matrix elementsll 
l^ts^fc/^cbl- Since and have been directly measured, one can combine 
these measurements to determine | Vts \ ■ In addition, the quantity Rk' provides 
information on the decay form factor in B — > K* We review in this section 
first the branching ratio B{B ^ +7) (and B{Hi, — > +7)) in the SM and 
then discuss estimates of |Vts| and Rk-- 



1.2 SM estimates of B[B ^ + 7) and B[Hb X,, + 7) 

The leading contribution to the decay 6 — > s + 7 arises at one-loop from the 
so-called penguin diagrams. With the help of the unitarity of the CKM matrix, 
the decay matrix element in the lowest order can be written as: 

M{b^s +j) = ^^Xt{F2{xt)~F2{xcWe''sa^.{mbR + m,L)& . (4) 

where Gf is the Fermi coupling constant, e — \/4,Traem, Xi — mf/rn^; i — 
u, c, t are the scaled quark mass ratios, and and are, respectively, -the 
photon four-momentum and polarization vector. The GIM mechanism El is 
manifest in this amplitude and the CKM-matrix element dependence is factor- 
ized in Af = VtbV^^. The (modified) InanrL-Lim fmiction ^2(2^1) derived from 
the (1-loop) penguin diagrams is given byQ: 

F2{x) = X {^x(?,x - 2) log a; ^ {x - \){%x^ + hx~ 7)] . (5) 

As the inclusive decay widths of the B hadrons are proportional to |Vc6p, 
the measurement of B(B — s- Xs 4- 7) can be readily interpreted in terms of 
the CKM-matrix element ratio At/|Vc6|. For a quantitative determination, 
however, QCD radiative and power corrections have to be computed, which 
we discuss next. 

The appropriate framework to incorporate QCD corrections is that of an 
effective theory obtained by integrating out the heavy degrees of freedom, 
which in the present context are the top quark and bosons. The effective 
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Hamiltonian depends on the underlying theory and for the SM one has (keeping 
operators up to dimension 6), 

Heffib s + 7) = -^A* J2 (6) 

where the operator basis, the lowest order coefficients Ci{mw) and the renor- 
malized coefficients Ci{^) can be seen elsewhere 0. The perturbative QCD 
corrections to the decay rate T{B — > Xs + 7) consist of two distinct parts: 

• Evaluation of the Wilson coefficients Ci(/i) at the scale /i = 0{mb). 

• Evaluation of the matrix elements of the operators Oi at the scale /i — 
Oimb). 

The Wilson coefficients are calculated with the help of the renormalization 
group equation whose solution requires the knowledge of the anomalous di- 
mension matrix in a given order in and the matching conditions, i.e., the 
Wilson coefficients Ci(/i — myy), calculated in the complete theory to the com- 
mensurate order. The leading logarithmic (LL) anomalous dimension matrix 
has been calculated by several independent groups First calculation of the 
next-to- leading order (NLO) anomalous-dimension matrix has been carried out 
by Chetyrkin, Misiak and MiinzE3. The matching conditions to order as have 
also been worked out in the meanwhile by several groups. Of these, the first 
six_corresponding to the four-quark operators have been derived by Buras et 
al.EJ, and the remaining two C7(/i = mw) and Cs{fi = mw) were worked out 
|by Adel and Yaollj. These lattetLiiave been recalculated by Greub and HiiEtiL 
113, confirming the earlier result 113. Recently, these matching conditions tHHl3 
have also been confirmed by Buras, Kwiatkowski and Potto. 

The NLO corrections to the matrix elements are of two kinds: 

• QCD Bremsstrahlung corrections b ^ sj + g, which are needed both to 
cancel the infrared divergences in the decay rate for B — > + 7 and in 
obtaining a non-trivial QCD contribution to the photon energy spectrum 
in the inclusive decay B ^ Xs + j- 

• Next-to-leading order virtual corrections to the matrix elements in the 
decay & ^ s + 7. 

The Bremsstrahlung corrections were calculated by Ali and Greub Ellll in the 
truncated basis (involving the operators Oi, O2, Ch-and Og) and jSubsequently 
in the complete operator basis by the same authorsll3 and by Pottlla. The NLO 
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virtual corrections were completed by Greub, Hurth and Wylerlil. These latter 
calculations have played a key role in reducing the scale-dependence of the LL 
inclusive decay width. All of these pieces have been combined to get the NT O 
decay width T{B ^ + 7) and the details are given in the literature E3t3. 

It is customary to express the branching ratio B{B ^ + 7) in terms of 
the semileptonic decay branching ratio B{B Xivi), 

B{B ^ Xsj) = [ ^iB^ l + Xs) y^ ^ ^^^^^^ 

i SL 



and the theoretical part can be expressed as 







A) . (8) 



Here, 



f{z) = l-8z + 8z-^-z*-12z'^lnz with z 



2 
2 

'^b,pole 



K(Z,fJ.b) \ mb^pole ) K(Z,fIb) V 3 TT / 

where rrib^poie i'mc,poie) is the b (c)-quark pole mass and rnb{^ — mi,) is the 
&-quark mass in the MS scheme. The function K(z,/2f,) represents the QCD 
correction to the semileptonic decay widthES, which depends on the scale jib = 
0{mb); its analytic fami is given by Nircll. The other functions {D and A) can 
be seen in literature £3. It should be remarked that while the Bremsstrahlung 
function A has been calculated in the complete operator basis, the virtual 
corrections contributing to D are still known in the truncated approximation. 
However, since the numerical values of the Wilson coefficients C3,...,Cq are 
very small, the contribution left out in the NLO expression in eq. (^) is expected 
to be small (not more than several percent). 

In addition to the perturbative QCD improvements discussed above, also 
the leading power corrections, which start in l/m^, have been calculated_t£ 
decay widths appearing in the numerator and denominator of eq. (H) c3c. 
The power corrections in the numerator have been obtained assuming that the 
decay B —>■ Xs+^ is dominated by the magnetic moment operator O7. Writing 
this correction in an obvious notation as 
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one obtains 6b = 1/2 Ai — 9/2 A2, where Ai and A2 are, respectively, the kinetic 
energy and magnetic moment parameters of the theacetical framework based on 
heavy quark expansion (henceforth caUed HQET)E3. Using Ai = —0.5 GeV^ 
and A2 = 0.12 GeV^, one gets Sb/ml ~ —4%. However, it turns out that 
the leading order (1/m^) power corrections in the heavy quark expansion are 
identical in the inclusive decay rates r{B Xg + 7) and r(i? — > Xiiyg), as 
far as Ai is concerned. The corrections proportional A2 differ only marginally. 
Thus, including or neglecting the l/m^ corrections makes a difference of only 
1% in the ratio (^) and hence in B{B —> + 7). 

Recently, the power corrections proportional to 1/m^, resulting from the 
intcp&pfpce of the operators O2 and O7 in B ^ Xg+'y, have also been worked 
outE30-E3. Expressing this symbolically as 

r{B^Xs + j) _^ , ^^^^ 



ro(B-.x, + 7) 
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one findsB SJml ~ +0.03. 

There exist several (marginally) differing numerical values of the branching 
ratio B{B ^ + 7) in the SM. Using \VtlVtb/Vcb\ = 0.976 ± 0.010 obtained 
from the unitarity constraints Ej, fixing the scale p,b = Wf,, varying the scale fib 
in the range 2mb > fJ^b > rnb/2, and using current values of the various input 
parameters, but ignoring the Sc/mj^ term, the short-distance contribution has 
been estimated by Chetyrkin et al. asli^: B{B Xs+7) = (3.28±0.33) x 10"^, 
where all the errors have been combined in quadrature. Including the Sc/rn?^ 
term, setting Jib — Mb ^-^id varying fib in the stated rangeJput keeping the othp; 
parameters the same as in the work of Chetryrkin et al. EJ, Greub and HurthE^ 
determine B{B ^ + 7) =X3.38 ± 0.33) x 10"''. Recently, a NLO value has 
been obtained by Buras et al. 0, yielding B{B X^+^f) = (3.48±0.3]±x 10"*. 
The shift in the central value from the one given by Greub and HurthEj is due 
to systematically discapdiiig the next-next-leading order terms, which were 
kept in the earlier workEjLJ, and the reduced error is due to treating the scale 
uncertainty in the numerator and denominator in eq. (||) independently. 

The numerical values for B{B — > ^37) quoted above Ej'oEj have been 
obtained by using the semileptonic branching ratio B{sl) = (10.4 ± 0.4)% 
taken from earlier data at T(4S'). This numbet,has changed somewhat in the 
meanwhile and the current measurements aretil: B{sl) — (10.49 ± 0.46)% (at 
T(4S')) versus B[sl) = (11.16 ± 0.20)% (at Z°). Updating B[sl) yields 

6(B^X, +7) = (3.51±0.32) X 10"* , (12) 

to be compared with the CLEG measurement B{B Xs-l-7) = (2.32±0.67) x 
10~^. The corresponding inclusive branching ratio at the Z'^ is obtained by 
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using in cq. the semileptonic branching ratio measured at the Z^. This 
gives 

B{Hb -^Xs+j)^ (3.76 ± 0.30) x lO"" , (13) 

to be compared with the ALEPH measurement B{Hi, — > + 7) = (3.29 ± 
0.98) X lO^**. The agreement between experiment and SM is good though the 
CLEO number is marginally (2(t) lower than the SM branching ratio. 

We would like to use the NLO SM-based theory and experiments to de- 
termine the CKM matrix element ratio |l^sVt6/146| and \Vts\- The two in- 
clusive measurements given in eqs. (^h and (|^), and the corresponding NLO 
SM-estimates given in eqs. ( p^ ) and p3| ) yield, 



-- 0.79 ± 0.11 (expt) ± 0.04 (th) (@T(4s)), 
ts"*": __ 0.91 ±0.14 (expt) ±0.04 (th) (@Z°). (14) 



Averaging the two measurements following the PDG prescription §1 gives the 
following weighted average for the CKM matrix element ratio: 



0.84 ± 0.09 (expt) ± 0.04 (th) 



=^ 0.84 ±0.10 (15) 

where the second row has been obtained by adding the theoretical and ex- 
perimental errors in quadrature. With the CKM unitarity, one has |— — 
|Vcs|; this equality holds numerically (within present precision) if one com- 
pares the l.h.s. obtained from the decay B —> Xg + 7 given in eq. (13) 
with the present determination of the r.h.s. from charmed hadron decayscE, 
lVcs\ = 1.01±0.18. Using the vahiaxif |V^fh| measured by the CDF collaboration 
§, I Vtbl = 0.99 ± 0.15 and notingB that \Vcb\ = 0.0393 ± 0.0028, finally yields 

l^tsl = 0.033 ±0.007 , (16) 

where all the errors have been added in quadrature. This is probably as direct 
a determination of \Vts\ as we will ever see, as the decay t ^ + s is too 
daunting to measure due to the low tagging efficiency of the s-quark jet. With 
improved measurement of B{B — > Xs + 7) and Vtb, one expects to reduce the 
present error on \ Vts\ by a factor of 2, possibly 3. 

The exclusive-to-inclusive ratio Rk-' has been worked out in a number of 
models. This involves estimates of the matrix elements of the electromagnetic 
penguin operator, implicitly assuming the SD-dominance. Taken on their face 
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value, different models give a rather large theoretical dispersion on Rk* . How- 
ever, one should stress that QCD sum rules, models based on quark-hadron 
duality, and the irrmroved lattice-QCD estimates of the 1997 vintage (being 
discussed by LynnE^) are theoretically more reliable. Concentrating only on 
them, some representativejesults are: 

Rk- = 0.20 ± 0.06 [Bah 0], Rk* = 0.17 ± 0.05 [Colangelo et al. 0] , 
Rk' = 0.16±0.05 [Ah, Braun & Simma S] , Rk- = 0.16±0.05 [Narison 0] , 
Rk' = 0.13 ± 0.03 [Ah & GreubEi] , Rk- = 0.16l^;J53 [Flynn El]. These 
estimates are consistent with each other and with the CLEO measurement 
Rk* = 0.181 ± 0.06. Summarizing this section, it is fair to conclude that SM 
gives a quantitative account of data in electromagnetic penguin decays in in- 
clusive rates, yielding a first determination of \Vts\ with an accuracy of ±20%. 
Further, both B{B ^ As + 7) and the ratio Rk* are in agreement with the 
dominance of the short-distance physics in these decays. 



2 Photon energy spectrum in _B As -I- 7 

Calculation of the photon energy spectrum is somewhat intractable as the 
main work-horse, namely HQET, does not quite make it to the very end of this 
spectrum. There is no alternative at present but to model the non-perturbative 
effects. We review the present state of the art. 

The two-body partonic process 6 — > S7 yields a photon energy spectrmn 
which is just a discrete line, l/{T)dr{b — > 57) = i5(l — x), where the scaled 
photon energy x is defined as a; = 2 E-ffri},/ {rn^ — m^). The physical photon 
energy spectrum is obtained by convoluting the non-perturbative effects and 
the perturbative QCD corrections, such as the ones arising from the decay 
h ^ s"f + g. The latter gives a characteristic Bremsstrahlung spectrum in x 
in the interval [0, 1] peaking near the end-points, — > i?™"^ (or a; — > 1) and 
— > (or x — *■ 0), arising from the soft-gluon and soft-photon configurations, 
respectively. Near the end-points, one has to improve the spectrum obtained 
in fixed order perturbation theory. This is done in the region a: — + 1 by 
isolating and exponentiating the leading behaviour in aemas(M)™ log^"(l ^ ^) 
with m < n, where /i is a typical momentum in the decay B — > As + 7. In 
this region, which is dominated by the magnetic moment operator O7, the 
spectrum can be symbolically expressed as 



' 3tt 



(17) 
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where C is a normalization constant (r = ms/mi,). 

C={1- rf{l + ,^ ^i(^--f^^iPole)\ ^,^^^^^, (^3) 

The running of ckg is a non-leading effect, but as it is characteristic pf 
QCD it modifies the Sudakov-improved end-point photon energy spectrum E3 
compared to its analogue in QEDEI. The expressions for il,i(x, r), ^2{x, r) and 
their deruaytiyes (denoted by primes) can be seen for non-zero s-quark mass in 
liter aturcHlaO. The expressions (with r — *■ 0) are: 

VLi = -21n^(l-.T)-71n(l-x), 

^2 = 10(l-.T) + (l-.T)2-|(l-a;)3-(l--a;)(3 + a;)ln(l-a;),(19) 



where the double logarithmic term in fii is universal La. The other terms are 
specific to the decay B ^ Xs + ^ ■ 

As long as the s-quark mass is non-zero, there is no collinear singularity 
in the spectrum. However, parts of the spectrum have large logarithms of the 
form as \og{rn'l / m^) , which are important near the end-point a; — > and are, 
in principle, present for any photon energy and should be resummed. The 
order Us corrected distributiniL for low photon energies is dominated by the 
operator Og and is given bjllII'Ej: 

U>o = <-^- a r(a;,r), (20) 



dx Stt 9 



where 



^, , /4-|-4r ^ ^\\-x + rx 

T{x, r) = 4 + 2a; In 

\x ~ rx J r 

(1 - x) [8 - (1 - r)a;(16 - 9x -|- 7rx) + {I - rYx^{l - 2x)\ 
xil — a; -|- rxy 

The difference between the fixed order spectmm and its exponentiated version 
can be seen in the work of Kapustin et al. Ej. If low energy photons can be 
detected in B ^ -t-7 (say, for <1 GeV), then such an experiment could 
help measure C^^^ . Unfortunately, with the SM values of the Wilson coeffi- 
cients, the partial branching ratio for B + 7 for low photon energies is 
too small to be measured even at B factories. However, with an anomaLpusly 
large C^'^^ , as has been entertained in the literature in other contextscOo, the 
C'l''^ ■^-dependent part of the photon energy spectrum may get appreciably en- 
hanced. It is worthwhile to measure the spectrum in the intermediate energies 
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(1.0 GeV < 2.0 GeV) to search for the effect of such anomalously enhanced 
Cg'^'^ contribution in beyond-the-SM scenarios. 

Implementation of non-perturbative effects is at present a model depen- 
dent enterprise and data ace, not, precise enough to distinguish various models 
proposed in the literaturell3'E3E3. We shall confine ourselves to the discussion 
of the photon energy spectrum calculated in a simple model, in which the b 
nuark in the B hadron is assumed to have a Gaussian distributed Fermi motion 
l3 determined by a non-perturbative parameter, pp. This model describes well 
the lepton energy spectrum in semileptonic decays B Xivi and it has also 
received some theoretical support in the HQET approach subsequently. 

The photon energy spectrum based on this model, including the QCD 
perturbative improvements, has been used both by the CLEOEl and ALEPH 
collaboration □ in the analysis of their data on B ^ Xs + 7.pAn analysis of 
the CLEO photon energy spectrum has also been undertaken IIZI to determine 
the non-perturbative parameters of this model, namely mi,{pole) and pp. The 
latter is related to the kinetic energy parameter Ai defined earlier in the HQET 
approach. The minimum of the CLEO data is obtained for pp — 450 
MeV and mi,{pole) ~ 4.77 GeV. However, the ilcr errors on these quantities 
are large (a sipular conclusion has been drawn in terms of Ai and mi,{pole) 
by Li and Yue3). The interesting question here is to determine if the non- 
perturbative aspects in the decays B Xlvt and S ^ can be described 
in terms of a universal shape function. To pursue this further requires lot more 
data which we hope will soon be forthcoming. 



3 Inclusive radiative decay B + 7 and constraints on the CKM 

parameters 

The quantity of interest in the decay B — > Xd + ^ is the high energy part of the 
photon energy spectrum, which has to be measured requiring that the hadronic 
system Xd recoiling against the photon does not contain strange hadrons to 
suppress the large- E'^ photons from the decay i? — > -(- 7, which is now 
the largest background. Assuming that such an experiment is feasible, one can 
determine from the ratio of the branching ratios B{B ^ Xd+^) /B{B Xg+'y) 
the parameters o£-the CKM matrix (in particular p and rj in the Wolfenstein 
parameterizationEJ) . 

In close analogy with the B Xs + 7 case discussed earlier, the complete 
set of dimension-6 operators relevant for the processes b d'j and b djg 
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can be written as: 

AG 



Heff (6 ^ d) ^ - -1 6 E (a^) (m) , 

.7 = 1 



(22) 



where — Vjb V*^ with j = u, c, t. The operators Oj, j — 1, 2, have implicit in 
them CKM factors. We shall use the Wolfenstein parametrizationc2l, in which 
case the matrix is determined in terms of the four parameters A, X ~ sin^c p 
and ?7, and one can express the above factors as : 

^u^AX^p~tr^), ec = -^A^ et^-en-Cc (23) 

We note that all three CKM-angle-dependent quantities £,j are of the same 
order of magnitude, 0{X^). It is ccpyenient to define the operators Oi and O2 
entering in Tief / (b — + d) as foUowsllj: 

01 = -^icLf3l'^bLa)idLaltJ.CLf3 - ^{uLf}^^bLa){dLalfi^ULp), 

02 = -^icLal''''bLa){dLpltJ.CLf3) " ^{uLal^bLa){dLI3l t^ULp) , (24) 

with the rest of the operators {Oj\ j = 3... 8) defined like their counterparts 
Oj in Heff{b — > s), with the obvious replacement s d. With this choice, the 
matching conditions Cj{mw) and the solutions of the RG equations yielding 
Cj (/i) become identical for the two operator bases Oj and Oj . The branching 
ratio B{B Xd + 7) in the SM can be generally written as: 

B{B^ Xd + -i) = DiX^ 

{{l-pf+r,^-{l-p)D2-7iD^ + Di}, (25) 

where the functions Di depend on various parameters such as mt,mb,mc, p,, 
andpqs. These functions were calculated in the LL approximation some time 
ago Lj and since then their estimates have been improved lB, making use of 
the NLO ca|k|Ulatipns discussed earlier. We shall assume, based on model 
calculations c§E3e3 and the 1 / power corrections discussed earlier in the 
context of B ^ Xg + 7, that the LD contributions are small also in B{B — > 

To get an estimate of B{B — > Xd + 7) at present, the CKM parameters p 
and r] have to be constrained from the unitarity fits, which yield the following 
ranges (at 95% C.L.)^: 

0.20 < 77 < 0.52, 
-0.35 < p< 0.35 . (26) 
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The (nominally) ppferred CKM-fit values at present are (p, 77) = (0.05,0.36), 
for which one getsE3 

6(B ^ Xd + 7) = (1.63±0.16) X 10-^ (27) 

where the error estimate follows from the one for B{B ^ + 7). Allowing 
the CKM parameters to vary over the entire allowed domain, one gets (at 95% 
C.L.) 

6.0 X IQ-^ < B{B Xd + -f)< 3.0 x 10"^ (28) 

The present theoretical uncertainty in this rate is a factor 5, which shows that 
even a modest measurement of B{B Xd + 7) will have a very significant 
impact on the CKM phenomenology. To the best of our knowledge, there is 
no experimental bound available on B{B ^ Xd + 7), but we hope that this 
decay will be measured in future at the B factories and CLEO. 

4 CKM-suppressed exclusive decays B{B ^ y + 7) 

Exclusive radiative B decays S ^ y + 7, with V = K* , p, uj, are also poten- 
tially very interesting for the CKM phenomenology EJ. Extraction of CKM 
parameters would, however, involve a trustworthy estimate of the SD- and 
LD-contributions in the decay amplitudes. We have argued that the decays 
J3 ^ + 7 and (B±,S") {K*'^,K*^°) + 7 are consistent with the dom- 
inance of the SD-contribution. There exist good reasons to believe that also 
the CKM-suppressed exclusive radiative decays are dominated by SD-physics, 
though one has to work out the LD-contribution on a case-by-case basis. More 
importantly, data on the various charged and neutral B meson radiative de- 
cays can be used directly to put meaningful bounds on the LD-contributions. 
Hence, despite skepticism in some quarters and in this conferenceEj, we believe 
that exclusive radiative B decays are worth measuring. 

The SD-contribution in the exclusive decays {B^,B^) {K*^,K*^) +7, 
{B^,B^) — > {p^,p^) + -f, B^ ^ UJ + ■y and the corresponding Bg decays, 
Bs —^ (f) + J, and Bg —> K*^ + 7, involve the magnetic moment operator Oj 
and the related one obtained by the obvious change s — > d, O7. The transition 
form factors governing the radiative B decays B ^ V + 'j can be generically 
defined as: 

{V,\\^i><j^,q''b\B) = ze^,p,e(^)p>^F#-^(0). (29) 

Here V is a vector meson with the polarization vector e*-^-' and stands for 
the field of a light u,d or s quark. In (p9|) the QCD renormalization of the 
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''P'^tivffb operator is implied. Keeping only the SD-contribution leads to obvi- 
ous relations among the exclusive decay rates, exemplified here by the decay 
rates for {B^,B°) ^ p + 7 and {B^,B°) ^ K* + j: 

r((i3±,S0)^(p±,p0)+7) 



r{{B±,B<^)-.{K*±,K*")+j)-''"''[\Vts\\ ' ^^"^ 

where Hi = [FsiBi p-f)/Fs{B, -> K*-/)]'^, which is unity in the SU{S) 
limit. (This is not being recommended— as the S'f7(3)-breaking effects have 
been calculated in a number of papers EjES.) Likewise, assuming dominance of 
SD physics gives relations among various decay rates 

r(s± ^ p±7) = 2 r(BO ^ PO7) = 2 r(BO ^ ^7) , (31) 

where the first equality holds due to the isospin invariance, and in the second 
SU{S) symmetry has been assumed. 

The LD-amplitudes in radiative B decays from the light quark intermediate 
states necessarily involve other CKM matrix elements. In the CKM-suppressed 
decays B V + j they are dominantly induced by the matrix elements of the 
four-Fermion operators Oi and 02- Estimates of these-cantributions have been 
obtained in the light-cone QCD sum rule approach US Using factorization, 
the LD-amplitude in the decay B^ + 7 can be written in terms of the 

form factors Fl" and F2, 

Along = ^^V^i,V:Jc2 + ^C^mpe^;'hi'>^ 



V2 ^wn-^^wy ".pc^ 

X { - ■ P) - pV] • (q^) + e^''"Vg/3 • 2^^2^((Z^)}32) 

The two form factors are found to be numerically close to each other in the 
QCD sum rule approach, F^ ~ = hence the ratio of the LD- and the 
SD- contributions reduces to a numberEa 

A, / A , — K>-B*^P*7 ^"faKTd 

■J^long I ■^short — J^l^l g ' -j,^ ■ \00 ) 



where 



± 



Rf,^" ^ = '^^ „ ' ' ■ — ^ = -0.30 ± 0.07 . 34 

The analogous LD-contributions to the neutral B decays pj and 
i?" — > uj^ are expected to be much smaller: 

L/s ^£^^_oi3^oo5 (35) 
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where the numbers are based on usingL^ 02/01 ~ 0.27±0.10 and ed/cu = —1/2 
is the ratio of the electric charges for the d- and u-quarks. This would then 

yield Rf"^''' ^ = 0-05- 

To get a ball-park estimate ofj-the ratio Along / Ashort, we take the central 
value from the CKM fits, yieldingca |yu6|/|Vt(j| — 0.33, which in turn gives. 



\'^long I ■^short\ — \-'^L/S I U/ ^r I ~ ' 

' Vtd Vtb 



A P"' 
long 

short 



< 0.02. (36) 



That the LD-effects remain small in 5° — > p-f decay has also been supported 
in an analysis basedjOp the soft-scattering of on-shell hadronic decay products 
j^o _^ pOpO _^ though this paper estimates them somewhat higher 

(between 4% and 8%). 

The relations (^l]), which obtain ignoring LD-contributions, get modified 
by including the LD-contributions to 

2r(B0 ^ p^) ~ 2r(i?o^a;7) " ^ ^ """^ ' ^ ^ 

where (Rl/s = -Rf/s"^ 



A{Rl/s) = 2 • RL/sVud _ ^-j2 + ^2 + (^L/s)'K?d _ + ^2 



(38) 



The ratio of the CKM-suppressed and CKM-allowed decay rates for charged 
B mesons gets likewise modified due to the LD contributions. Following earlier 
discussion, we ignore the LD-contributions in T{B K*^). The ratio of the 
decay rates in question can therefore be written as: 



\'[il-pf+fj']{l + AiRL/s)) (39) 



The effect of the LD-contributions is modest but not negligible, introducing an 
uncertainty comparable to the ^ 15% uncertainty in the overall normalization 
due to the S'C/(3)-breaking effects in the quantity 

Neutral B-meson radiative decays are less-prone to the LD-effects, as ar- 
gued above, and hence one expects that to a good approximation (say, better 
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than 10%) the ratio of the decay rates for neuti-al B meson obtained in the 
approximation of SD-dominance remains vaHd LJ: 



^KdX'[{i-pr+fi'], (40) 



where this relation holds for each of the two decay modes separately. This 
ratio is at par with the mass-difference ratio AM^/AMg in the neutral B- 
meson sector, as both involve a reliable estimate of the 5'[/(3)-breaking effects 
but otherwise reflect the dominance of the SD- physics. i— in 
Finally, combining the estimates for the LD- and SD-form factors e3E3, 
and restricting the Wolfenstein parameters in the allowed range given earlier, 
yields 

S(B±->p±7) = (1.5±1.1) X lO-*' , 

= (0.65 ± 0.35) X IQ-^ , (41) 

where we have used the experimental value for the branching ratio B{B — > 
K* + 7) El. The large range reflects to a large extent the poor knowledge of the 
CKM matrix elements and hence experimental measurements of these branch- 
ing ratios will contribute greatly to determine the Wolfenstein parameter p and 
77. Present experimental limits (at 90% C.L.) areB: B{B^ P^l) < 1.1x10"^, 
S(B" p7) < 3.9 X lO"'"^ and B{B° ^7) < 1.3 x 10"^ The constraints on 
the parameters (p, ry) following from them are, however, not yet competitive 
to the one foUowingJ'rom unitarity and lower bound on the mass difference in 
the B° - 5° sector E. 

5 Inclusive rare decays B Xs£^i^ in the SM 

The decays B Xsi^i~ , with £ =-e^jj.,,T., nmvide new avenues to search for 
physics beyond the standard model Ea'Ca^EllHEj. The branching ratio B{B — > 
Xs + 7) constrains the magnitude of Cj'^ but the sign of C^'^ is not deter- 
mined by the measurement of B(B — s- Xs + 7). This sign is in general model 
dependent. It is known that in SUSY models, both the negative and positive 
signs are allowed as one scans over the allowed SUSY parameter space. The 
B Xst^t~ amplitude in the standard model (as well as in several exten- 
sions of it such as SUSY) depends on the coefficient C^^ and additionally on 
the coefficients of two four- Fermi operators, Cg and Ciq. It has been arguedE3 
that the signs and magnitudes of all three coefficients C^^ , Cg and Cio can, in 
principle, be determined from the decays B ^ Xg + ^ and B Xsi'^£~ by 
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measuring the dilepton invariant mass and Forward-Backward charged lepton 
asymmetry E3. 

The SM-based rates for the decay b si'^l" , calculated in the_free quark 
decay approximation, have been known in the LO approximation Ed for some 
time. The NLCLcontribution reduces the scheme-dependence of the LO result 
in these decaysE3. In addition, long-distance (LD) effects, which are expected 
to be very important in the decay B Xsi^i^ , have also been estimated from 
data on the assumption that they arise dominantly due to the charmonium res- 
onances J/ip^ijj'^ ■■■ through the decay chains B XsJ/tp{ip', ■■■) Xs£^£~ ■ 
The resulting dilepton distribution is then a coherent sum of the resonating 
(LD) and mildly varying (SD) contributions. Data near the resonances can 
be used to better parametrize the LD contribution in future than is the case 
now. Recently, these LD-contributions have also been predicated upon assum- 
ing that far from the resonant-region they canJpe determined in terms of the 
1/m^ contributions in the HQET approachESH. This remains an interesting 
conjecture but impossible to test experimentally, as, first of all, they repre- 
sent just a class of power corrections and, more importantly, this theoretical 
framework breaks down near the resonances. 

The leading (1/mf,'^) power corrections to the partonic decay rate and 
the dileptop invariant mass distribution have been calculated in the HQET 
approachcJ; these resjilts have, however, not been confirmed in a recent inde- 
pendent calculation^, which finds that the power corrections in the branching 
ratio B{B Xsi'^£~) are small (typically —1.5%). The corrections in the 
dilepton mass spectrum and the FB asymmetry are also small over a good part 
of this spectrum. However, the end-point dilepton invariant mass spectrum is 
not calculable in the heavy quark expansion and will have to be modeled. As 
an alternative, non-perturbative effects in S ^ Xs£~^£~ haszE also been esti- 
mated O, using the Fermi motion model discussed earlier a. These effects, 
which model power corrections in l/mt, are also found to be small over most 
of the phase space except for the end-point dilepton mass spectrum where they 
change the underlying parton model distributions significantly and have to be 
taken into account in the analysis of dataEil. 

Taking into account the spread in the values of the input parameters, 
/i. A, TOt, and Bsl discussed in the previous section in the context of B{B 
s -|- 7), the following branching ratios for the SD-piece have been estimated 

B{B^Xse+e-) = (8.4 ± 2.3) x 10-^ 
BiB ^ X,fi+^i-) = (5.7± 1.2) X 10-^ 
BiB~^XsT+T-) = (2.6±0.5) X 10"^ (42) 
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where theoretical errors and the error on Bsl have been added in quadrature. 
The experimental upper limit for the inclusive branching ratio forJiie decay 
B was quoted by the UAl collaboration some time agoL3, B{B 

Xs^jL^ llZ) < 5.0 X 10^^. This limit has been put to question in a recent CLEO 
paperB. From the CLEO data, a hmit B{B < 5.8 x IQ-'^ has 

been set on this decay, with B{B XgC+e^) < 5.Z^ 10^^. Combining the di- 
electron and di-muon upper hmits, CLEO quotesEJ B{B Xs£^(^) < 4.2 x 
10~^ (all limits are at 90% C.L.). This is a factor 6 away from the SM estimates. 
As far as we know, there is no experimental limit on the mode X^r+r". For 
a more detailed discussion of the modes B Xs£'^£^ , some related exclusive 
decays such as Bj-.^^ {K, K*)£~^£~ , and other rare i?-decay modes, we refer 
to recent reviews eIeZI. The CKM-suppressed decays B Xd£~^i~, which are 
expected to be typically a factor 20 below their corresponding CKM-allowed 
decay rates, and thcip-jrolc in determining the CKM parameters have been 
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